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We elucidate the physical mechanisms by which electromagnetic excitations in planar chiral waveg- 
uides exhibit broken time reversal and parity symmetries as recently observed in metallic nanos- 
tructures (Phys. Rev. Lett. 91, 247404 (2003)). Furthermore, we show that concomitant with 
these broken symmetries the electromagnetic excitations can acquire fractional winding numbers, 
revealing analogies between light scattering from planar nanostructures and anyon matter. 
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A recent experimental study of optical interactions 
with non-magnetic metal planar chiral nanostructures 
consisting of arrays of four-fold chiral gammadions pro- 
vides intriguing evidence of broken time reversal (T) 
symmetry Planar chiral objects are characterized 

by the fact that they cannot be brought into congruence 
with their enantiomeric or mirror image form without 
lifting the object from the plane. Among the intrigu- 
ing properties of chiral arrays is that they display non- 
reciprocal polarization effects in diffraction which resem- 
ble the well-known non-reciprocity of the Faraday effect 
0. In the experiment of Ref. chiral nanostructures 
whose topographies were mirror images of each other 
were found to lose there mirror symmetry when viewed in 
polarized light, an effect that also occurred for individual 
chiral gammadions. It has been proposed that these un- 
usual symmetry properties result from the chiral optical 
response associated with plasmon and volume modes in 
the grooves of the structure Q . 

In this paper we theoretically show that electromag- 
netic (EM) excitations in planar chiral nanostructures 
can exhibit broken T and parity (V) symmetries, and a 
physical argument is given for how these broken symme- 
tries can be manifested in the properties of light scattered 
from the nanostructures. We also point to analogies be- 
tween light scattering from planar chiral nanostructures 
and from anyon matter . 

To describe the experiment of Ref. 0] we consider a 
thin metallic (gold-titanium) film with an array of four- 
fold chiral gammadions carved in it. We neglect coupling 
between individual gammadions based on the high losses 
suffered by propagating surface plasmon-polariton waves 
propagating between gammadions due to the titanium 
layer covering the structure. These losses, however, are 
negligible for EM excitations confined in the grooves of 
individual gammadions as considered here, since the ti- 
tanium layer is absent in the grooves. We thus consider 
an individual four-fold chiral gammadion centered on the 
origin, an example of which is shown in Fig. ^ Here the 
gammadion lies in the (x, y) plane and we have chosen 




FIG. 1: [colors online] An AFM image of a fragment of the 
planar chiral structure showing one gammadion carved in a 
thin gold/titanium film, as described in Ref. i2j. The elec- 
tromagnetic excitations considered in the present paper are 
concentrated in the grooves of the structure that form a single 
gammadion element. The characteristic sizes of the gamma- 
dion's elements, their width and depth are 1.5 fim, 700 nm, 
and 120 nm, respectively. 



to align the x-axis with the central portion of one gam- 
madion groove and the y-axis with another, the z-axis 
being perpendicular to the plane of the gammadion. For 
future use we also introduce a cylindrical coordinate sys- 
tem for which points in the plane of the gammadion are 
represented by the radial coordinate p and azimuthal an- 
gle (p as indicated in Fig. ^ Then, for an EM field 
of frequency lo the dielectric tensor around the gamma- 
dion structure has elements e^i/(r), with p,,i> = x,y,z, 
where we omit the explicit lo dependence for simplicity 
in notation. We approximate the tensor as Hermitian 
and its diagonal elements as real by the virtue of the 
fact that for the optical frequencies of interest in gold 
Rsi^fifi) 1^1 IiT^i^fifj.) \, and the four-fold gammadion 
symmetry implies €xx — ^yy The generally complex off- 
diagonal tensor element e^y — ty^ describes polarization 
coupling in the (x, y) plane due to the chiral gamma- 
dion patterning, and we assume that this is the domi- 
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nant source of polarization coupling and hereafter set all 
other off-diagonal tensor elements to zero. To proceed 
we expand the tensor elements exyi'^ into azimuthal ba- 
sis functions 



(1) 



The basis functions cxp{i£(j)) are useful for characterizing 
the chirality of the structure as each has an associated in- 
teger winding number £ such that the phase varies from 
zero to 2'Ki upon circling around the z-axis. The real- 
ity of the diagonal dielectric tensor elements implies that 
they have equal weighting of the winding numbers ±£, 
£^2 = meaning that the chirality is not mani- 

fest in the diagonal elements. In contrast, the generally 
complex off-diagonal tensor elements can manifest the 
chirality of the gammadion. For illustrative purposes we 
consider the simple case that the off-diagonal elements 
are dominated by a single winding number i, the contri- 
bution from —£ being set to zero, for example £ ^ 4 for 
a gammadion with four-fold symmetry and £ — —4 for 
the corresponding enantiomeric form, and we write the 
off-diagonal elements as exy{r) — {^yx{r))* — k{p, z)e'^'^. 
Finally, for the diagonal elements we ignore, to lowest- 
order, the patterning of the gammadion, so that the di- 
agonal elements ^^^{z) are functions of z only, so the 
dielectric tensor becomes 



^xx{z) K{p,z)e^^'^ 
e(f) = I Kip,z)e-''^ exx{z) 
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The dielectric tensor Q describes essentially a non-local 
medium for which the off-diagonal tensor elements are 
space-dependent through the chiral terms exp(±i^0), and 
we have taken the function z) to be real so as to 
highlight the effects of the chirality. We further assume 
that |e^^| >> 

To proceed we examine EM excitations that are bound 
to an individual chiral gammadion or waveguide. These 
EM waves are time-harmonic solutions of the Maxwell 
equations in the presence of the gammadion waveguide 
with electric fields of the form E{r, t) = [E{r^ oj, t)e~*"* -I- 
c.c.J/2. In the slowly-varying envelope approximation 
when the vector electric field envelope E{r,uj,t) changes 
slowly on the time-scale 1 /uj and we neglect second-order 
time derivatives, Maxwell's equations lead to 



-V X V X + 2iw/ioe(f) 



dE 
'dt 



poioh{r)Eir)^0, (3) 



with the generally complex dielectric tensor e(r) given 
by Eq. (0). For transverse electric (TE) waves whose 
polarization lies in the (x, y) plane we find, using the 
dielectric tensor (jSJ and neglecting the off-diagonal ele- 
ments to dominant order, that V • D = V • {e{r)E) = 



ixx{z)\7E = 0, which yields V • ^ = 0. This restric- 
tion to TE polarization is justified if we assume that the 
diagonal dielectric tensor element £xx{z) provides suf- 
ficient confinement along the z-axis so that the usual 
designation of TE polarization (electric field perpendic- 
ular to z) and transverse magnetic (TM) polarization 
(magnetic field perpendicular to z) for waveguides ap- 
plies For strong confinement along the z-axis we 
assume that the field structure along the z-direction is 
dominated by the normalized fundamental mode u{z) 
determined from [d?/dz'^ -\- fiouj'^exx 
with /3" the corresponding eigenvalue. Then substi- 
tuting E{r,uj,t) = u{z)[ex£x{p,z,t) + ey£y{p,z,t)] in 
Eq. ©, with Cx^y denoting unit vectors along the re- 
spective directions, projecting out the mode structure 
along z, and introducing the spinor notation £{p, <j), t) = 
col {£x{p,4',t),£y{p,(l),t)), Eq. Q reduces to the follow- 
ing effective planar equations for the EM excitations cou- 
pled to the chiral waveguide 



where the Hermitian Hamiltonian is given by 
pQLU^R{p)e''^'^ 



H = Ht 



poLu^R{p)e-''^''' 



(4) 
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with diagonal part i/j^ = — (V^^ + / being the unit 
(2 X 2) matrix and V'^ the transverse Laplacian acting 
in the (x^y) plane, and k{p) = J c?z|u(z)pK(p, z). We 
remark that by projecting out the mode structure u{z) 
along z we have rendered the system effectively two- 
dimensional (2D). However, the physical system as de- 
scribed by Eqs. Q is still sensitive to the selection of 
the direction of the z-axis since the gammadion winding 
number £ is defined with respect to the angular direction 
of phase advance around the z-axis of the chiral term 
exp («£(/)). The key ingredient for the following discussion 
is that we are considering TE waves for which the elec- 
tric field has only two vector components instead of three, 
and this allows for our two component spinor approach. 
This 2D planar approximation for TE waves is valid as 
long as the structure is tightly confined and single-mode 
along the z-axis. We note that Eqs. are analogous to 
the Bogoliubov equations for electrons in a superconduc- 
tor 6], with k(p, z) playing the role of the pair potential. 

We next examine the time-reversal properties of the 
spinor wave Eq. (0J. For our system a suitable antilinear 
time-reversal operator is T = KqI, where Kq is the com- 
plex conjugation operator . By operating T on Eq. Q 
we obtain the time-reversed equation 
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{T£) =THT~\T£), 



and for H given in Eq. (jSJ we find 

THT~^ =H* = VHV^^. 



(6) 



(7) 
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For our 2D reduced system described by the spinor Eq. 
Q the hnear parity operator V has the effect of reversing 
(j) —(/), which can be implemented for the Hamiltonian 
using V = cr^, with cr^ being the PauH-spin matri- 
ces. Since a structure and its corresponding enantiomeric 
form are interconverted by space inversion Eq. O 
tells us that time-reversal is equivalent to spatial inver- 
sion. This feature was noted in Ref. 0] and simply means 
that under time-reversal incident EM waves approach the 
planar structure from the opposite direction and perceive 
it as the enantiomeric structure of opposite chirality. In 
the absence of chirality £ = 0, when the Hamiltonian 
in Eq. is purely real, the time-reversed and space 
inverted Hamiltonians in Eq. ((JJ) are both equal to H 
meaning that both T and V symmetries are intact. In 
contrast, in the presence of chirality T and V symmetries 
are simultaneously broken for the planar chiral structure. 
However, Eq. imphes that {VT)H{VT)-^ = H, so 
that the system is always invariant under the combined 
action of VT symmetry jj]. 

More generally, the off-diagonal matrix element e^j,(r) 
of the dielectric tensor may expressed as a sum over wind- 
ing numbers £ as in Eq. , a special example of which 
is the single dominant winding number case considered 
above. Then, unless eiy = (eiy^^)* for all i and e(r) is 
real, both T and V symmetries are simultaneously bro- 
ken. Essentially the time reversal and parity asymmetries 
in planar chiral structures are underpinned by the non- 
locality of the gammadion response, meaning that exci- 
tation of the extended gammadion as a whole is crucial. 
A planar structure consisting of individually radiating 
atoms could not produce such a symmetry breaking in- 
teraction. 

We are now in a position to discuss how these new re- 
sults relate to the experiment of Ref. which provided 
strong experimental evidence of broken T and V symme- 
tries. Distinct polarization phenomena for media where 
T and V symmetries are simultaneously broken have pre- 
viously been discussed in relation to anyon matter in 
high-temperature superconductors, initially by Wen and 
Zee and later by Canright and Rojo 0. For the pur- 
pose of analyzing how the broken symmetries applicable 
to TE excitations may lead to observable effects, it is 
important to realize that an incident EM field can cou- 
ple energy into the planar TE waves of individual chiral 
waveguides as strict momentum conservation is relaxed 
at the edges of the grooves of the gammadions, and by the 
same argument EM energy can be out-coupled from in- 
dividual chiral waveguides. Furthermore, previously the 
in-plane patterning of the diagonal dielectric tensor el- 
ement exxiz) was neglected. At next order the radial 
dependence of ea;a;(/9, 0, z) provides a grating structure 
that can in-couple and out-couple energy to radiation 
modes analogous to the emission mechanism of surface 
emitting circular grating lasers |l5l |. This coupling be- 
tween the incident, out-coupled, and TE waves of the 



chiral waveguides, alongside the losses present in real sys- 
tems, give the TE waves a finite lifetime. Furthermore, 
these mechanisms offer an explanation of why the fields 
scattered from chiral planar nanostructures can manifest 
the broken T and V symmetries applicable to the bound 
TE waves, since the diffracted fields from the nanostruc- 
tures have their origin in energy out-coupled from the TE 
waves excited by the incident field. Thus, it is physically 
reasonable that the scattered EM fields will inherit the 
broken T and V symmetries exhibited by the planar TE 
waves. 

We have also found that concomitantly with the break- 
ing of the T and V symmetries the corresponding TE 
excitations bound to chiral waveguides can acquire frac- 
tional winding numbers and we now demonstrate this 
interesting property which also elucidates the vector na- 
ture of the TE waves. Inspection shows that consistent 
stationary solutions of Eq. Q), corresponding to time- 
harmonic solutions at frequency cj, assume the form 

where the winding number of the x-polarized component 
of the field is (m -I- £) , and m is that for the y-polarized 
component, m = 0, ±1, ±2, . . . being an integer that we 
use to label the solutions. This form of solution ensures 
consistency between the azimuthal variation of the elec- 
tric field polarization components in the presence of the 
phase dependent coupling terms exp(±j^(/)) in the dielec- 
tric tensor. We introduce the average winding number 
for a TE vector field as 

< Wem >- ■^0^'' ^'^ P'^P^^ ^)^^^^^ (P' ^) 

(9) 

where we have defined the EM winding number operator 
Wem = —i{d/d(f>)I. For the EM wave solutions in Eq. 
© we find 

<Wem >=m + e{l + si)/2, (10) 

where 

Iq'' d(t>f^ pdpllaxip)]"^ -\ay{p)\^] 

Sl = — 2^ — . (tij 

lo'' d(l>J^ pdp[\ax{p)\^ + \ay{p)\^] 

is the spatially averaged and normalized value of the 
element si = — l^^y!^) of the Stokes vector s = 

{si, 82,83) 01 over the plane of the individual waveg- 
uide. In general, — 1 < si < 1, so that the EM winding 
number in Eq. IjlUII can be non-integer, that is, frac- 
tional. The average winding number of the field is as- 
sociated with its orbital angular momentum that has its 
origin in the spatial structure of the fields 0|, as dis- 
tinct from the spin angular momentum associated with 
the polarization state of the field. The average spin angu- 
lar momentum is proportional to the spatially averaged 
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value of the S3 = 2Im{£x{p,4>,t)*£y{p,(j),t)) component 
of the Stokes vector i which for £ 7^ yields S3 = for 
the EM solutions in Eq. (|SJ) by virtue of the azimuthal 
integrals over (j). This implies that in the presence of pla- 
nar chirality the state of polarization of the vector field 
components varies over the plane in such a way that the 
average spin angular momentum is zero, and the aver- 
age angular momentum of the coupled fields is equal to 
the orbital angular momentum associated with the frac- 
tional modal winding number. To see the consequence 
of this we realize that here we are discussing the classi- 
cal TE vector modes of the structure, and we anticipate 
that upon quantization of the field we can add photons 
to each of these vector modes . Associated with each 
photon added to vector TE mode of index m there will 
be an average angular momentum directed along the z- 
axis < Jz >= h < Wem >= h{m + e{l + si)/2), that is, 
the average angular momentum of the quantized vector 
modes can be fractional. The average winding number 
for the EM excitation can be fractional since the devia- 
tion from the integer values expected for a bosonic field is 
fully accounted for by the angular momentum associated 
with the chiral waveguide J^] . 

It is worthwhile commenting on the nature of the 
EM mode interaction with the chiral gammadion. As 
noted earlier, polarization is associated with the spin of 
the EM field, whereas the spatial variation of the EM 
field describes orbital angular momentum, for example, 
a Laguerre-Gaussian laser field [13. From Eq. © we 
see that for £ the relative phase between the two 
vector field components varies with azimuthal angle (j), 
meaning that the polarization state of the field is not 
spatially homogeneous, as is obvious from the polarized 
images provided in Ref. flj. Furthermore, since the two 
vector field components have different winding numbers 
they carry different orbital angular momenta. Thus, the 
chiral waveguide intertwines the spin and orbital angu- 
lar momenta for the EM modes in such a way that they 
cannot be separated, in contrast to optical activity, for 
example, in which the molecular chirality causes polar- 
ization rotation but otherwise leaves the orbital angular 
momentum unchanged. 

In summary, we have elucidated the physical mecha- 
nisms by which EM waves coupled to chiral waveguides 
can display broken T and V symmetries. Symmetry is 
recovered in the form of enantiomeric reversibility VT 
based on the equality of excitations in the forward and 



time-reversed processes involving left- and right-handed 
structures, respectively. We have argued that these bro- 
ken symmetries should also be manifest in the fields scat- 
tered from the structures as reported in Ref. 0. Fur- 
thermore, we have shown that the TE excitations bound 
to chiral waveguides can acquire fractional winding num- 
bers. We finish by drawing attention to the analogy be- 
tween the light scattering from anyon matter and chiral 
nanostructures, since both involve broken T and V sym- 
metries and excitations with fractional angular momen- 
tum 
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